Minimal size of a barchan dune 
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Barchans are dunes of high mobility which have a crescent shape and propagate under conditions 
of unidirectional wind. However, sand dunes only appear above a critical size, which scales with 
the saturation distance of the sand flux [P. Hersen, S. Douady, and B. Andreotti, Phys. Rev. Lett. 
89, 264301 (2002); B. Andreotti, P. Claudin, and S. Douady, Eur. Phys. J. B 28, 321 (2002); G. 
Sauermann, K. Kroy, and H. J. Herrmann, Phys. Rev. E 64, 31305 (2001)]. It has been suggested 
by P. Hersen, S. Douady, and B. Andreotti, Phys. Rev. Lett. 89, 264301 (2002) that this flux 
fetch distance is itself constant. Indeed, this could not explain the proto size of barchan dunes, 
which often occur in coastal areas of high litoral drift, and the scale of dunes on Mars. In the 
present work, we show from three dimensional calculations of sand transport that the size and the 
shape of the minimal barchan dune depend on the wind friction speed and the sand flux on the 
area between dunes in a field. Our results explain the common appearance of barchans a few tens 
of centimeter high which are observed along coasts. Furthermore, we find that the rate at which 
grains enter saltation on Mars is one order of magnitude higher than on Earth, and is relevant to 
correctly obtain the minimal dune size on Mars. 

PACS numbers: 45.70.-n, 45.70.Qj 



I. INTRODUCTION 

Dunes are beautiful sand patterns formed by the wind, 
and are found in deserts, and also along coasts. Dunes 
have different sizes and shapes which depend on the con- 
ditions of wind and sand. The simplest and best under- 
stood dunes are the barchans P, 0, H, 0, H, H, 0, H> H 
[Tol . [Ill [r2| . These dunes appear under conditions of uni- 
directional wind and low sand availability. They have a 
windward side, two horns, and a slip face at the lee side, 
where avalanches take place (fig. QJi) . Barchans are sub- 
ject of scientific and also environmental interest because 
of their high rate of motion. For instance, barchan dunes 
1 — 5 m high may cover 30 — 100 m in a year. Such dunes 
are found very often on coasts, where they emerge from 
the sea sand, the grains of which, after being deposited 
onto the beach, dry and are thereafter transported by 
the wind. Yet, barchans nucleate only when a sand heap 
reaches a minimal size,a process which is still poorly 
understood [H [H El (H . 

The grains that form sand dunes are carried by the 
wind through saltation [l[: When the shear velocity u* 
of the wind exceeds the threshold friction speed for sand 
transport, grains are lifted out of the sand bed and en- 
trained into motion. These grains are next accelerated 
downwind, impacting after a certain distance I back onto 
the ground (fig. [Tjp). After a grain-bed collision, other 
grains may be ejected in a splash [lfl [lT], [l8| from the sur- 
face, which yields a cascade process in which the number 
of grains in saltation increases exponentially. However, 
the momentum transfer from the air to the grains yields 
a deceleration of the wind ("feedback effect" [19() which 
leads to flux saturation (fig. [Tb). When the slope of the 
surface exceeds the angle of repose, 9 r « 34°, a slip face 




FIG. 1: a. Sketch of a barchan dune showing the windward 
side, horns and slip face. In the inset we see the definitions 
of dune width W and length L. b. Schematic diagram with 
the main elements of saltation. The number of ejected grains 
(indicated by dashed arrows) is proportional to the velocity 
of the impacting grain, W m P [ lal. which is of the order of the 
average grain velocity (v) [141 ]. The mean saltation length of 
the grains is defined as £. c. Behaviour of the sand flux q as 
a function of horizontal distance x. The characteristic length 
of flux saturation is called saturation length, £ s . 



is developed. 

The characteristic distance to reach flux saturation is 
called saturation length £ s (fig. [lb), and determines the 
minimal size of dunes. The saturation length increases 
with the length of the saltating trajectories, £, which in 
turn scales with ^drag = ^Pgrain/pfluid, where d is the 
mean diameter of the grains, /C gr ain is their density, /9fl u id 
is the density of the driving fluid, and ^drag is the distance 
within which a sand grain lifted from the bed reaches the 
wind velocity [l3| . The proportionality between the min- 
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FIG. 2: Barchan dunes in Morocco. 
Wide Web. 
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imal dune size and the saltation length and ^drag has been 
verified in the field and also in the laboratory. Hersen et 
al. (2002) 10], supported by previously reported mea- 
surements on several terrestrial barchan dune fields, ob- 
served that the smallest dunes have heights of the order 
of 

^dragj and that the minimal width Wmin of a barchan 
dune is approximately 20 times £drag- Since the sand of 
aeolian dunes is constituted by quartz grains of mean di- 
ameter d = 250 fim. ctncl density Pgrain 

= 2650 kg/m 3 , 

we obtain, using the air density pfl U id = 1.225 kg/m 3 , 
^drag (and therefore H min ) ~ 50 cm, and W mm around 
10 m (fig. It is somewhat intriguing that Hersen et 
al. (2002) could find a similar relation for the min- 
imal size of aquatic dunes although saltation in water is 
much attenuated because of the high density of the driv- 
ing fluid [l(|. In their experiment, they found dunes 800 
times smaller than the aeolian dunes, since pfl u id of water 
is 10 3 kg/m 3 . 

As observed recently, there are other physical variables 
besides £ri rag which appear relevant to explain the scale of 



dunes [15j. Barchans are also found on Mars, where the 
atmospheric density is 100 times smaller than on earth. 
Because the mean grain diameter d = 500 /im [2(| and 
density p gra i n ~ 3200 kg/m 3 are not very different from 
the earth's, we would expect dunes on Mars to be ac- 
cordingly 200 times larger, and at least around 2 km in 
width, since £drag on Mars is nearly 100 m. However, we 
find on Mars dunes of only a few hundred meters in width 
(fig. [3]). On the basis of this surprising inconsistence, it 
has been suggested that still unknown microscopic prop- 
erties of the martian saltation must be understood and 
taken into account in order to correctly predict the size 
of dunes on Mars [15j . 

Recently, much understanding has been gained 
through dune modellization, which encompasses the main 
ingredients of the physics of dunes and accounts for flux 
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FIG. 3: Mars Global Surveyor (MGS) Mars Orbiter Cam- 
era (MOC) images of Mars dunes, a. Barchan dunes in the 
Arkhangelsky Crater (42°S, 25°W). In b., c. and d. we see 
north polar dunes in (76.3°N, 263.5°W), (76.6°N, 255.9°W) 
and (71.7°N, 51. 3° W), respectively. We see that very small 
dunes do not display slip face or horns. 



saturation [1J,|21|. This model reproduces quantitatively 
the shape of barchan dunes measured in the field, and 
correctly predicts the sand flux and wind profile over 
dunes, their velocity and also the scale dependence of 
their shape [13, [23j]. In the present work, we use the 
dune model to calculate how small a barchan can be, 
and we find that the size and the shape of the minimal 
dune depend on the wind shear velocity and on the sand 
flux at the interdune 0, [HIM 111 H ■ Furthermore, we 
present the first three dimensional calculation of dunes 
on Mars. 

The dune model and the calculation procedure are de- 
scribed, respectively, in sections II and III. In section IV 
we present and discuss our results. In Section V, we ex- 
tend the calculations to Mars. A summary of our main 
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conclusions is given in section VI. 



II. THE DUNE MODEL 

Here we give a brief presentation of the dune model, 
and refer to refs. [3, HH, [23| for a more extensive de- 
scription. 

The fundamental idea of the model is to consider the 
bed-load as a thin fluid-like granular layer on top of an 
immobile sand bed. This picture is based on the original 
idea of Bouchaud et al. (1994) (2?| who proposed a model 
for the motion of grains during avalanches on sand piles, 
which has been latter adapted for the modelling of dunes 
[2l| . The dune model combines an analytical description 
of the average turbulent wind velocity field above the 
dune with a continuum saltation model which allows for 
saturation transients in the sand flux. 



A. Wind shear stress 

Sand transport takes place near the surface, in the tur- 
bulent boundary layer of the atmosphere [2a ]. In this 
turbulent layer, the wind velocity u(z) at a height z may 
be written as: 



u* z 
u(z) = — In — , 

K Z 



(1) 



where k = 0.4 is the von Karman constant, it* is the wind 
shear velocity, which is used to define the shear stress r = 
PfluidU^, and zo is the aerodynamic roughness, u* and zq 
are two independent variables which may be determined 
experimentally. One way to obtain them is to measure 
the wind velocity at different heights z, and to plot the 
data as a linar-log curve u(z) vs logz. The inclination of 
the straight line of the fit is the shear velocity it* , and 
the value of z for which u(z) — is the roughness zq. 
This method has been applied for instance to determine 
the wind profile, shear velocity and surface roughness at 
the Pathfinder landing site on Mars [2q . 

For sand transport to occur, the wind shear velocity 
must exceed a threshold value u*t which is, on Earth, 
around 0.22 m/s (25|. This minimal shear velocity de- 
pends on the grain diameter, the density of the grains and 
of the fluid, gravity, and also on the packing of grains [28| . 
and is predicted to be almost ten times larger on Mars 
[29j . On the other hand, the aerodynamic roughness zq is 
still subject of research [3(| [M|. It is distinguished from 
the roughness Zg and which is of the order of a few tens of 
microns, and which is due to the microscopic fluctuations 
of the sand bed when the grains are at rest, zq means 
the "apparent" roughness which is a consequence of the 
motion of saltating grains. Bagnold (1941) [l| already ob- 
served that zq must be larger than ZQ and , and increases if 
there are pebbles or rocks. A value of zo close to 1.0 mm 
has been often reported from measurements of saltation 



on a sand bed [25(, while on Mars zq is larger, around 1.0 
cm [H]. 

According to eq. |T]), the wind velocity over a flat 
surface increases logarithmically with the height above 
the ground. A dune or a smooth hill can be considered as 
a perturbation of the surface that causes a perturbation 
of the air flow over the hill. In the dune model, the 
shear stress perturbation is calculated in two dimensional 
Fourier space using the algorithm of Weng et al. [32j 
for the components t x and r y , which are, respectively, 
the components parallel and perpendicular to the wind 
direction. The following expressions hold for the Fourier- 
transformed shear stress perturbation components: 



7~x (kx i ky) 



2 h(k Xj ky)h^ 



\k\uHi) 

2 ln(£|fc s |) + 4r + 1 + isign(^)7r 



and 



2 h{k x , ky)k x ky 



(2) 



(3) 



where the axis x(y) points parallel (perpendicular) to 
the wind direction, k x and k y are wave numbers, \k\ = 

kl + fc2, T = 0.577216 (Euler's constant) and C is the 

characteristic length of a hill [331 ] . It is defined as the 
horizontal distance between the crest, which is the posi- 
tion of maximum height -ff max , and the position of the 
windward side where the height is i/ max /2. The variable 
C is computed iteratively, i.e. it is not a constant param- 
eter but depends on the size of the hill at each iteration. 
U{1) = u(l)/u(h m ) is the undisturbed logarithmic profile 
(fTJ) calculated at height I, which is given by 



2k 2 C 
In l/zo ' 



(4) 



and normalized by the velocity at the reference height 
h m = C/^/logC/zq, which separates the middle and up- 
per flow layers [33(. The shear stress in the direction i 
(i = x,y) is then given by: 



n = i [t (i + n 



(•5) 



where tq is the undisturbed air shear stress over the flat 
ground. From the shear stress, the sand flux is calculated 
according to the continuum saltation model (l4| . 

In what follows, we give a brief presentation of the 
sand transport equations and refer to Sauermann et al. 
(2001) [H[ and Schwammle and Herrmann (2005) [13] 
for the extensive derivation of the saltation model. 



B. Continuum saltation model for sand transport 

The saltation model is derived from the mass and mo- 
mentum conservation in presence of erosion and external 
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forces. The sand bed represents an open system which 
can exchange grains with the saltation layer, and the ero- 
sion rate T at any position (x, y) on the bed must balance 
the local change of the sand flux, V-<f 

The erosion rate is defined by the difference between 
the vertical flux of grains leaving the bed and the rate <f> 
at which grains impact onto the bed: 



r = 0(n-l), 



(6) 



is defined as the length 

: v e J e (2(v)/g), where v e J e 



The mean saltation length £ 
of a ballistic trajectory [14| : £ = 
is the initial vertical velocity, g is gravity and (v) is the 
average grain velocity, which is independent of the shear 
velocity u„, and is calculated as explained below. Fur- 
ther, v c z iC — aAwhor, where a is an effective restitution 
coefficient for the grain-bed interaction [f4T | . In this man- 
ner, £ — (l/r)[2(v) a/g], where r = (v)/Avh m -- In this 
manner, the saturation length is then written: 



where n is the average number of splashed grains. The 
flux of saltating grains reduces the air born shear stress r a 
("feedback effect"). At saturation, the number of ejecta 
nearly compensates the number of impacting grains (n = 
1), and r a is just large enough to sustain saltation, i.e. 
r a is close to the threshold r t = Pf\ u au1 t [Hj]. In this 
manner, we write n as a function ?i(T a /T t ) with n(l) = I. 
Expansion of n into a Taylor series up to the first order 
term at the threshold yields 



IT) 



where 7 = dn/d(r a /r t ), the entrainment rate of grains 
into saltation, determines how fast the system reaches 
saturation [l4| . The parameter 7 depends on microscopic 
quantities of the grain-bed and wind-grains interactions, 
which are not available within the scope of the model. 
Therefore, 7 must be determined from comparison with 
measurements or microscopic simulations. 

The rate <fi at which the grains impact onto the bed is 
defined as |<f|/4 where I is the average saltation length. 
Substituting this expression for ip and eq. ([7]) into eq. 
(J6j> , the balance r = V-cf yields 



\q\ 



But r a = t — Tg, where r g is the contribution of the grains 



to the total shear stress 
the above expression: 



141. In this manner we rewrite 
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The grain born shear stress is defined as r g = 0Ai>h O r, 
where Awh or gives the gain in horizontal velocity of the 
particle after one saltation trajectory. Therefore, r g = 
A?;hor|<7|/A and the above expression for the sand flux 
may be written as 



r-t [ Av hoi /£ 

1 - m 

T-Tt 



where we can identify two important physical quantities: 
the saturation length £ s = [^/ 7 ]T t /(r — r t ), and the satu- 
rated sand flux q s = (r — Tt)£/ Avh or - 
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while the following expression is obtained for the satu- 
rated sand flux q s : 



2a(v) 
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(9) 



The resulting equation for the sand flux is a differential 
equation that contains the saturated flux q s at the steady 
state (eq. |H]) and the saturation length £ s (eq. [5]): 



v. 9 -=i|,i(i-M 

4 V <is 



(10) 



The mean velocity of the saltating grains, (v), is deter- 
mined from the balance between the drag force acting on 
the grains, the loss of momentum when they impact onto 
the ground, and the downhill force. To calculate (v) we 
need indeed to take into account the modification of the 
air flow due to the presence of the saltating grains. How- 
ever, the model equations do not account for the complex 
velocity distribution within the saltation layer. Instead, 
a reference height Z\ is taken, between the ground at the 
roughness height ZQ and = d/20 [34| and the mean salta- 
tion height z m , at which the "effective" wind velocity in 
the equilibrium, u ff, is calculated [34] ]: 



u* t I , Zj 
Ueff = < In — — T + 2 




-1-1 




Next, the grain velocity, v, is calculated numerically from 
the equation [Til ]: 



3 Pfluid Cd 

4 Pgr a in d 



(u eS - v)\Ueff ~ v\ 



gv 
2a\v\ 



gVh = 0. (12) 



The grain velocity obtained from eq. (fT2j) is in fact the 
average grain velocity at the steady state, i.e. \v\ = (v), 
since u e ff is the reduced wind velocity after flux satu- 
ration has been achieved fli| . In the simple case of 
the two-dimensional flow over a sand bed, we obtain 
(v) = u c ff — Vf I y/2a, where Vf is the grain settling velocity 
[34| . The drag coefficient Cd and the model parameters 
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a, z m and z\ are obtained from the expressions derived 
inref. H. 

On the other hand, we can obtain the constant from 
measurements of the saturation length. Sauermann et al. 
(2001) have determined it from reported measurements 
and microscopic simulations of saltation [H, [H, H3] , and 
found r7 « 0.2 [l4j|. So far we could not estimate r or 7 
separately and it is the quantity r"f which can be obtained 
from comparison with experimental data. However, a 
simple argument shows that, while 7 increases with the 
amount of splashed grains, r should not be much different 
on Mars [3f|. 

The steady state is assumed to be reached instanta- 
neously, since it corresponds to a time scale several orders 
of magnitude smaller than the time scale of the surface 
evolution. Thus, time dependent terms are neglected. 
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FIG. 4: The initial surface is a gaussian hill of height H and 
characteristic length a. 



exceeds 34°, then the slip face is formed and the separa- 
tion streamlines are introduced. Calculations consist of 
the iterative computation of steps (i) — (iv). 

III. CALCULATIONS 



C. Surface evolution 

The surface is eroded wherever the sand flux increases 
in the direction of wind flow, and sand deposition takes 
place if the flux decreases. The time evolution of the 
topography h(x,y,t) is given by the mass conservation 
equation: 



dh 
dt 



1 



(13) 



where p S and = 0.62p gra j n is the mean density of the im- 
mobile sand which constitutes the sand bed If 
sand deposition leads to slopes that locally exceed the 
angle of repose 6* r » 34°, the unstable surface relaxes 
through avalanches in the direction of the steepest de- 
scent. Avalanches are assumed to be instantaneous since 
their time scale is negligible in comparison with the time 
scale of the dune motion. 

For a dune with slip face, flow separation occurs at 
the brink, which represents a discontinuity of the sur- 
face. The flow is divided into two parts by streamlines 
connecting the brink with the ground. These stream- 
lines define the separation bubble, inside which eddies 
occur [2lj. In the model, the dune is divided into slices 
parallel to wind direction, and for each slice, the separa- 
tion streamline is identified. Each streamline is fitted by 
a third order polynomial connecting the brink with the 
ground at the reattachment point as described by Kroy 
et al. (2002) [2l|. Inside the separation bubble, the wind 
shear stress and sand flux are set to zero. 

The simulation steps may be summarized as follows: 
(i) the shear stress over the surface is calculated using the 
algorithm of Weng et al. (1991) [32]; (ii) from the shear 
stress, the sand flux is calculated using eq. (fTUj). where 
the saturation length £ s and the saturated sand flux q s 
are calculated from expressions (0) and ([9]), respectively; 
(iii) the change in the surface height is computed from 
mass conservation (eq. ([131)) using the calculated sand 
flux; and (iv) avalanches occur wherever the inclination 



We perform calculations with open boundaries, and 
uni-directional wind of friction speed and a constant 
sand influx q- m at the inlet. The influx is interpreted 
as the average interdune flux in a dune field. The dune 
size depends on the initial volume of sand, but its shape 
does not depend on the initial topography. We use the 
simplest initial surface, which is a gaussian hill of height 
H and characteristic length a as shown in fig. [4] The 
gaussian shape evolves until it achieves a final barchan 
shape, which is a consequence of the shear velocity 
and the influx q- m . 

The shear velocity it* may have v ery different values 
depending on the location of the field 25, 39]. Illenberger 
and Rust (1988) found that rates of sand deposition can 
be ten times higher in coastal dune fields than those in 
mid-desert sand seas, due to the abundant sand supply 
on sandy beaches and the higher energy of coastal winds 
[4(i| . In our calculations, the important quantity is the 
shear velocity associated with the typical value of u*/u*t 
conditioned by u*/u*t > 1.0. The typical u*, is thus, 
an average over the values of wind friction speed above 
the threshold u*t- While the typical u* is around 0.32 
m/s {u*/u*t ~ 1.45) at the barchan field in Qatar (6j, it 
is close to 0.39 m/s (u*/u*% w 1.8) at Jericoacoara and 
Lengois Maranhenses in northeastern Brazil [22], |4^] • 
On the other hand, the interdune flux has been subject 
of field measurements S 24 1 and may vary significantly 



even in the same area [24|. It depends on the amount 



of loose sand available between dunes, on the size and 
distribution of large immobile particles, humidity and the 
presence of sparse vegetation. The flux between barchans 
in a field is normally smaller than 50% of the maximum 
flux [13 . 

Here we investigate for the first time using three di- 
mensional calculations of barchan dunes how the field 
variables influence the minimal dune size. In particular, 
we look for the smallest barchan dune for each set of 
{it*/u*t, 9in/?s}- Below this minimal dune size, the hill 
has neither a slip face nor horns, and is called a dome. 
We also want to investigate how the shape of the mini- 
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mal dune depends on the wind strength and the flux in 
the intcrdunc area. In our calculations, we set a constant 
value of sand influx at the inlet, which is a small frac- 
tion of q s and is considered, for simplicity, homogeneous 
along the y axis (perpendicular to sand transport). To 
determine the minimal size, we perform simulations of 
dunes of different sizes, where we begin with sand hills 
of different heights H with a constant H/a « 0.2. The 
dimensions of the largest dune below which slip face and 
horns are absent define the minimal dune size. 

This procedure finds immediate application in Plane- 
tary Science. Figure [3] shows images of Mars dunes sent 
by the Mars Global Surveyor (MGS) Mars Orbiter Cam- 
era (MOC). We see that Mars barchan dunes have differ- 
ent shapes, and that the size and the shape of the domes 
and of the dunes are particular for each field. In fact, the 
properties of Mars dunes are related to the specific condi- 
tions of wind and flux in the field [ll[ . Although circula- 
tion models have been applied to estimate the strength of 
martian winds, very few data from in situ measurements 
are available, and no estimation has been reported for 
areas with barchans on Mars. It is also not possible to 
determine the sand flux which forms Mars dunes since 
rovers never landed on barchan dune fields. Modelliza- 
tion of Mars dunes may help significantly to understand 
the wind behaviour on Mars, and may also yield informa- 
tions about microscopic properties of the martian sand 
transport. 

In the next section, we present the results of our calcu- 
lations for terrestrial barchan dunes. We consider quartz 
grains with diameter d — 250 /im, since this is a repre- 
sentative value for the sand of barchan dunes 25]. The 
other quantities which we need to solve the model equa- 
tions are the angle of repose of the sand, # r = 34°, 
and the density of the grains, p gl . a i n = 2650 kg/m 3 ; the 
gravity, g = 9.81 m/s; the density of the driving fluid, 
Pfiuid = 1-225 kg/m 3 ; the dynamic viscosity of the air, 
n = 1.8 kg/m-s, which is used to calculate the trajectories 
of the grains [341 ] : and the minimal wind shear velocity 
for saltation, u*t = 0.22 m/s. 



IV. RESULTS AND DISCUSSION 

The shear velocity m* is the only of the studied field 
variables which explicitly enters expression ([8]) for the 
saturation length of the flux, £ B . We expect the dimen- 
sions of the smallest dune obtained in calculations to de- 
crease with the shear velocity, since the saturation length 
also decreases with it*. We calculate the minimal width 
Wmin and length L m i n of the barchans for different values 
of between u*% and 2.3u*t- We notice that a value of 
2.3 u*t ~ 0.5 m/s is associated with a wind velocity of 
8.6 m/s or 31 km/h at a height of 1 m, using a roughness 
length Zq = 1 mm. Wind velocities larger than 8.0 m/s 
are among the strongest measured on dune fields, and 
are found for instance in northeastern Brazil [4l|, |42j . 

As we can see in the main plot of fig. 03 the minimal 
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FIG. 5: Main plot: Minimal dune width Wmm, as a function 
of the relative shear velocity M*/ii„t, obtained with gi n /g s = 
0.20. The lower inset shows the "excentricity" I/ m i n /W m i n . 
In the upper inset, we show L m i n (squares) and W m in (dia- 
monds) as functions of the characteristic length of flux satu- 
ration, £ s , calculated with the corresponding values of M*/tt*t 
in the main plot. The straight lines are displayed to show the 
linear increase of the dimensions of the minimal dune with 
the saturation length £ s . 



dune width decreases from 13.0 m to 3.0 m if it* changes 
from 1.37u*t or 0.3 m/s to 2.3it*t or 0.5 m/s. On the 
basis of the result of fig. 03 it is interesting that the 
smallest barchans found in the Lengois Maranhenses in 
northeastern Brazil [52] have widths between 5 and 10 
m. Wind shear velocities reported for that region reach 
values between 0.35 and 0.45 m/s (1.6 and 1.9 times it*t) 

[HEllil- 

The result of the main plot of fig. [5] is a consequence 
of the expression obtained for the saturation length (eq. 
([S]) ), which is proportional to the mean saltation length £ 
but also decreases with the pre-factor 7[(u*/it* t ) 2 — 1] of 
the multiplication process of saltation. The denominator 
of eq. ([8]) gives the efficiency of the wind in carrying 
grains into saltation. This increases, in turn, with the 
relative wind strength (r — T t )/r t and with the amount 
of grains available from grain-bed collisions — the faster 
the population of saltating grains in the air increases, the 
faster the flux saturates, and the shorter is the saturation 
length £ s . The linear increase of the minimal dune width 
Wmin and length L m j n with £ s can be seen in the upper 
inset of fig. El where we see that the saturation length 
is the relevant length scale of barchan dunes. The width 
of the smallest barchan, VF m i n , is around 12 — 14 times 
£ Sl while the minimal dune length L nnn is approximately 
22-24 4- 

In the lower inset of fig. [5] we see that the shape of 
the dome also changes with the shear velocity. We call 
^min/Wmin the "excentricity" of the smallest dune. Fur- 
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FIG. 6: In the box on top, we show dunes of width 6.5, 12, 21 
and 48 m, from the left to the right, calculated with u* = 0.36 
m/s and qin/qs = 0.2. Threshold shear velocity is u*t = 0.22 
m/s. The main plot shows the linear relation between L and 
W for different values of shear velocity: u, = 0.3 (circles), 
0.36 (squares), 0.41 (triangles) and 0.46 m/s (stars). The 
upper-left-hand-corner inset shows the excentricity L/W as 
a function of the dune width W for the corresponding values 
of «*. In the bottom-right-hand-corner inset, we see that the 
width W of the dune increases with the height H in the same 
manner for all values of u* . 



ther, we also see in fig. [5] that dunes of different sizes 
have different excentricities L/W. Therefore, the quan- 
tity Amin/Wmin is particular to the minimal dune, and is 
useful for the characterization of a dune field. 

Therefore, we see that the strength of the wind plays a 
relevant role for the minimum size of dunes. Here it must 
be emphasized that the wind strength we consider is a 
mean value of the wind velocities (above the threshold) 
associated with sand transport. In a given field, the wind 
is normally very fluctuating in time, a large fraction of 
which its velocity may be even below the threshold for 
saltation (2f| . Time series of wind speed measured in the 
field may be found for instance in refs. [43|,[44j]. The wind 
velocity u* we use in our calculations is interpreted as a 
representative value of winds strong enough for saltation 
to occur — an average over the wind velocities above the 
threshold. 

Indeed, even in the same field, dunes may display dif- 
ferent shapes. Dunes with asymmetries, for example, 
may be a result of small, local fluctuations in the wind 
direction and/or variations in the topography jTl . Fur- 
ther, the shape of the dune depends on the amount of 



sand influx, as we call the flux at the interdune area for 
the downwind barchan, which may vary significantly in 
a dune field [p. 

One factor which determines the sand influx is the dis- 
tribution of dunes upwind. While the net flux just after 
the slip face of a barchan dune is zero, the flux is nearly 
saturated at the tip of dune horns. The sand which leaves 
the horns is transported through the interdune area and 
reaches the windward side of the downwind dune. The 
amount of sand transported depends on several variables. 
First, it depends on the shape of the dune horns; dunes 
with thicker horns lose more sand. We notice that this 
property can only be captured in a three-dimensional cal- 
culation of dunes, since the horns are not included in a 
two-dimensional model [H, Hl|. On the other hand, the 
sand flux which arrives at the dune depends on the trans- 
port properties of the interdune area. While interdune 
areas with vegetation and humidity may trap the sand, 
on bedrock, the flux is essentially conserved, while the 
sand flux on areas with much sand is saturated. 

Here we investigate the effect of the amount of incom- 
ing sand on the dune shape and the minimal dune size. 
First, we show in fig. [7^ how different a barchan dune 
of width 180 m appears with different values of shear 
velocity and sand influx q m /q s from 1 to 50%. Kroy et 
al. (2005) UH have already shown that the aspect ratio 
H/L increases with u*. Furthermore, there is an interest- 
ing feature which we can only see with three-dimensional 
calculations. The "slim" shape of barchans is character- 
istic for areas of low q- lri /q Sl and as the influx increases, 
the dunes become "fat" . The differences between "slim" 
and "fat" barchan shapes have been noticed by Long and 
Sharp (1964) g[, who reported that these dunes also be- 
have differently. Long and Sharp (1964) observed that 
the "fat" barchans of the Imperial Valley, California, are 
the more "morphologically complex and areally larger 
dune masses" , and may be also a result of dune inter- 
action, coalescence and dune fusion. Their interpreta- 
tions are in accordance with our calculation results for 
a larger interdune flux and therefore a higher amount of 
sand (influx) at the upwind of barchans. 

In fig. [TJd, we show the shape of the smallest dunes 
obtained with different values of u*/it*t and q\ n /q s - Fur- 
thermore, fig. [5]shows the excentricity imin/Wmin of the 
smallest dune as a function of qi n /q s for different values 
of u*/u*%. In the inset of this figure, we show the mini- 
mal dune width Wmin as a function of q- m /q s . We see that 
the minimal dune width increases with the sand influx, 
and that the excentricity L m i n /W m - ln decreases almost 
linearly with q m /q s - In summary, we conclude that the 
shape of the smallest dune provides a local indicative of 
the amount of mobile interdune sand in a barchan field. 



V. MINIMAL DUNE SIZE ON MARS 

We have applied the dune model to study dunes on 
Mars. In the context of barchan dunes study, Mars 
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FIG. 7: a Barchan dunes of width W = 180 m calculated 
with different values of relative wind friction speed u*/u*t 
and influx qm/q s - The gray area in (b) has dimensions 73 x 
104 m, and shows the minimal dunes corresponding to each 
set {u»/u*t, <?m/q s }. 



constitutes itself a particularly interesting scene where 
winds are of much more uni-directional character than 
on Earth 45]. We see in fig. [3] that Mars barchans 
have many different shapes, which may indicate particu- 
lar conditions of wind and flux. This has been previously 
speculated by Bourke et al. (2004) [ll|, who performed 
a careful study of the morphology of dunes in craters 
and on the north pole from images of Mars dune fields. 
Bourke et al. (2004) measured the height, length and 
width of dunes and concluded that while on the north 
pole dunes have a clear correlation between these quanti- 
ties, intra-crater barchans display heterogeneities, which 
have been attributed to "topographically induced, vari- 
able wind regimes" and also to limited sand supply in 
craters [ll|. More generally, Bourke et al. (2004) ob- 
served that intra-crater dunes are longer and wider than 
the north polar dunes [III • 

The dune model equations make use of sand and at- 
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FIG. 8: Main plot: Excentricity L m i B /W m m of the smallest 
dune as a function of qi n /q s calculated for different values of 
shear velocity: u*/u*t = 1-37 (squares), 1.64 (circles) and 
2.50 (triangles). The corresponding values of minimal dune 
width VKmin are shown in the inset. 



mospheric properties as well as of field quantities of wind 
and flux, which may be estimated through comparison 
of the different dune shapes with calculation results . 
Therefore, using the dune model to study Mars dunes 
might yield a helpful tool in the understanding of the 
aeolian conditions and geologic history of Mars. Indeed, 
the model depends on the local temperature, which is 
found to vary between 150 and 300 K, and the atmo- 
spheric pressure, which varies between 5 and 10 mb on 
Mars. These values are necessary to calculate the den- 
sity and the viscosity of the air, which is used to obtain 
the trajectories of the grains [34(. The local conditions 
of the atmosphere have been repeatedly measured by the 
MGS Radio Science Team |47|, whose data are used in 
the calculations of the present work. 

However, there has been no previous estimation of the 
entrainment rate of grains into saltation on Mars, which 
we need to calculate the saturation length £ s . The rate 
of sand entrainment, 7, is proportional to the number of 
grains ejected after grain-bed collisions because saltating 
particles come mainly from the splash. Further, Ander- 
son and Haff (1988) [l6| have shown that the number of 
grains which are splashed from the surface scales with 
the velocity of the impacting grain, D; mp (fig. [T}d). In 
order to extend this relation to Mars, we must rescale 
this velocity with the escape velocity, \fgd, which is the 
velocity necessary to escape from the potential trapping 
at the sand bed surface [48J. This velocity is around 4.5 
cm s _1 both for Mars and for Earth. On the other hand, 
the impact velocity t>i mp scales with the average velocity 
of the saltating grains, (v) [Oj], and therefore we obtain 
7 cx (v)/y/gE. 

We remark that the grain velocity used in the scaling 
relation for 7 is the velocity in the equilibrium: for wind 
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TABLE I: Average velocity (v) of saltating grains on Earth 
and on Mars as a function of the relative shear velocity u*/u, t - 



velocities which are not much larger than the thresh- 
old both on Mars and on Earth — typical values of 
used in our calculations are a few units of M*t — we use 
the approximation that the ratio between the velocities 
of the grains on both planets, VMars/^Earth is roughly 

(«Mars)/ (wEarth)- 

Typical values of the average velocity of saltating 
grains on Mars and on Earth are shown in Table 1 as 
a function of the relative wind friction speed u»/u»t- We 
see that the grain velocity on Mars is one order of mag- 
nitude larger than on earth. Moreover, we see that the 
dependence of (v) on it, for both cases is less noticeable, 
which is a surprising result and is a consequence of the 
model equation (fT2]) , from which the grain velocity is cal- 
culated using the reduced wind velocity (eq. pip ) due to 
the presence of the saltating grains . From the scaling 
of the grain velocity, it follows that 7 is constant which 
is proportional to u*t, and is independent of u*. 

On the basis of the difference in grain velocities on 
Mars and on Earth, we obtain that 7 on Mars is one order 
of magnitude higher than on Earth. Thus, the constant 
rj which appears in the denominator of eq. ([5]) and is 
around 0.2 for saltation on earth [14| . has value i=s 2.0 for 
Mars. 

What is the consequence of a 10 times larger entrain- 
ment rate on Mars? The higher rate at which grains 
enter saltation on Mars amplifies the "feedback effect" 
and reduces the distance of flux saturation. Therefore, 7 
plays an essential role for the scale of dunes, in addition 
to the mean hopping length and to the wind velocity 
and sand flux as shown in the previous section. 

Summarizing, the following are the most important pa- 
rameters in the problem of modelling of Mars dunes: The 
atmospheric temperature T and pressure P; the grain 
diameter d and density /O gra in; and the gravity g. From 
T and P we determine the air viscosity and the den- 
sity, which are used to obtain the model parameters of 
saltation as mentioned above. Finally, we also obtained 
an estimation for the rate of entrainmcnt of grains into 
saltation on Mars, which has been the unknown micro- 
scopic quantity of saltation relevant for the calculation 
of the saturation length £ s . We next look for the values 
of it»/u*t and q-m/q s which reproduce the shape and the 
size of Mars dunes. 



250 
200 
S 150 
^ S 100 
50 




2.2 



12.1 



J 8 1.9 



0.1 



0.2 

q. Iq 

* in J c 



0.3 



1.5 



2.0 

u lu . 



2.5 



3.0 



FIG. 9: Main plot: Minimal width W m i n of barchan dunes on 
Mars as a function of the relative shear velocity u*/u*t. The 
value VK m i n ~ 200 m is associated with a shear velocity around 
1.45 u*t or 2.8 m/s. In the inset, we show the excentricity of 
the minimal dune as a function of qi n /q s , calculated using 
u, = 2.8 m/s. We see that the excentricity Z/min/Wmin ~ 2.0 
of the domes in the Arkhangelsky Crater is reproduced with 
an interdune flux of 20% of the saturated flux. 



We study the shape of the barchan dunes in the 
Arkhangelsky Crater on Mars. We have chosen these 
dunes because they have a wide spectrum of sizes, in- 
cluding two domes of width 200 m and length ss 400 
m which indicate the minimal dune size (fig. [3^,). The 
Arkhangelsky Crater is situated in the Southern Hemi- 
sphere of Mars, in a region where the mean atmospheric 
pressure is around 6 mb and the temperature is nearly 
200 K. These values are used in the calculations, and we 
also use the values of grain diameter and density listed in 
Section I, as well as the threshold shear velocity u*t = 2.0 
m/s, and the entrainment rate 7 ten times larger as ob- 
tained above. 

In the main plot of fig. [H] we show the minimal dune 
width TF m i n as a function of the relative shear velocity 
u*/u*t obtained with parameters for the Arkhangelsky 
Crater. In this figure, we show an astonishing finding: 
the minimal size of Mars dunes decreases from 250 m 
for it* close to the threshold shear velocity to nearly 20 
m for u* ss 6.0 m/s or three times the threshold. This 
is a variation of one decade in the minimal dune size. 
Furthermore, we conclude that the value of W m [ n = 200 
m is obtained with around 1.45 u+t or 2.8 m/s. Next, 
we show in the inset the excentricity of the minimal dune, 
^min/Wmin, as a function of qi n /q s . We see that the ratio 
^min/Wmin « 2.0 is obtained if the average interdune 
flux gi„ is approximately 20% of the saturated flux q s . 

Figure [10] shows the results obtained using u* = 
1.45 u* t and qi n /q s — 0.20. In fig. [TUa . we show three 
Arkhangelsky dunes of different sizes next to dunes cal- 
culated with the model, while in fig. [TUb we show the 
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FIG. 10: a. From the top to the bottom, images of Arkhangel- 
sky dunes of width ~ 200 m (domes), 450 m and 650 m, next 
to calculations of dunes of similar sizes, b. Length L versus 
width W of the Arkhangelsky barchans (circles) and of dunes 
calculated using u* = 2.8 m/s and qi n /q s = 20% (full line). 



length L as function of width W of the Arkhangelsky 
dunes (circles) and of the dunes obtained in calculations 
(full line). We see that the values of u*/u*t and qi n /q s 
obtained for the Arkhangelsky Crater on Mars not only 
reproduce the minimal dune but also describe well the 
dependence of the shape on the dune size. We also con- 
clude that the elongated shape of intra-crater barchan 
dunes on Mars .IjJ is associated to wind friction speeds 
close to the threshold shear velocity for saltation. 

We note that if the value of the entrainment rate 7 
were the same on Mars as on Earth, the scale of the 
dunes would be 10 times larger than observed [50J; i.e. if 
we neglected the influence of the larger splash on Mars, 
on the basis of fig. [Tp, in the computation of the satura- 
tion length. Kroy et al. (2005) reported the failure of the 
pure scaling £ s oc £dmg in a recent work where they stud- 
ied the scale of dunes on Mars [HI]. In their work, they 
found that this scaling overestimates the dune scale on 
Mars — as mentioned in Section I of the present work — 
and concluded that differences in sand transport mecha- 
nisms on Mars and on Earth could be associated with this 
discrepancy. Here we have shown that the larger splash 
on Mars is the missing link to understand the scale of 
Mars dunes. Furthermore, the role of the wind velocity 
for the dune size in other dune fields on Mars is subject 
of ongoing research [38j . 

The saturation length obtained using parameters for 
Mars is approximately 17 m, and the minimal dune width 
is nearly 12 times the saturation length. Moreover, the 
value of m* found for the Arkhangelsky Crater is within 
the range of estimated wind friction speeds, between 2.0 
and 4.0 m/s, at the Mars Exploration Rover Meridiani 
Planum [5l[ and Viking 1 [52j landing sites on Mars. The 
results presented here provide indeed indirect, theoretical 
evidence that the scarce atmosphere of Mars could still 
transport sand and form dunes in the present. 

VI. CONCLUSIONS 

We presented an extensive study of the minimal dune 
size using a set of three-dimensional equations of sand 



transport. Our calculations show that the shape and the 
size of the minimal dune depend on the wind velocity 
and on the local sand flux in the interdune area. We 
have shown that the minimal dune size decreases with 
the wind friction speed, and the excentricity (length over 
width) of the smallest dune is determined by the inter- 
dune sand flux for a given wind velocity. On the basis 
of our findings, we can explain the appearence of small 
dunes on the coasts where the average friction speed of 
winds above the threshold reaches values between 0.35 
and 0.45 m/s. 

Furthermore, we have found that the rate at which 
grains enter saltation on the Planet Mars is one order of 
magnitude higher than on earth, and has a strong impli- 
cation on the minimal size of Mars dunes. The higher 
entrainment rate amplifies the feedback effect associated 
with the loss of momentum of the wind due to the ac- 
celeration of the grains, and shortens the distance of flux 
saturation. As a consequence, the minimal size of Mars 
dunes estimated from the scaling with a constant flux 
fetch distance is reduced by a factor of 10. It would be 
interesting to use this new insight to make a full micro- 
scopic simulation for the saltation mechanism of Mars 
similar to the one that was recently achieved by Almeida 
et al. (2006) With the entrainment rate calculated 
for Mars, we could obtain the observed minimal dune size 
in the Arkhangelsky Crater. We could also give an esti- 
mate for the sand flux in the interdune area and found a 
wind friction speed of 2.8 m/s. This value is well within 
the range of u* previously estimated for Mars [5l|, l52j |. 
which likely correspond to maximum values of wind fric- 
tion speed that occur during dust storms. The range of 
wind velocity on the present Mars extends in fact to much 
lower values [53j, [5J] , which explains the lack of observa- 
tions of dune activity in the last martian decades. We 
can say that the value of u*/u*t = 1.45 found for Mars 
was the typical value of relative wind velocity at the time 
when the dunes were active. The model equations repro- 
duced moreover the dependence of the dune shape on the 
size using the variables of wind and flux obtained from 
the minimal dune size on Mars. 
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